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Abstract 

We propose a simple model of the neutrino mass matrix which can explain the solar 
and atmospheric neutrino problems in a 3(z/£+z/r) framework. Assuming that only two 
right-handed neutrinos are heavy and a Dirac mass matrix has a special texture, we 
construct a model with four light neutrinos. The favorable structure of flavor mixings and 
mass eigenvalues required by those neutrino deficits is realized as a result of the seesaw 
mechanism. Bi-maximal mixing structure might be obtainable in this scheme. Since it 
contains a light sterile neutrino, it has a chance to explain the LSND result successfully. 
We consider an embedding of this scenario for the neutrino mass matrix into the SU(5) 
grand unification scheme using the Froggatt-Nielsen mechanism based on U ( 1 ) ^ xU(l)^ 2 . 
Both a small mixing angle MSW solution and a large mixing angle MSW solution are 
obtained for the solar neutrino problem depending on the charged lepton mass matrix. 
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1. Introduction 

Recently the existence of non-trivial lepton mixing has been strongly suggested through 
the atmospheric and solar neutrino observations whose results can be explained by as- 
suming the neutrino oscillations |1|, 0, |3| . The predicted flavor mixing is much bigger than 
the one of quark sector. The explanation of this feature is a challenging issue for the con- 
struction of a satisfactory grand unified theory (GUT) and a lot of works have been done 
0, In most of them the smallness of the neutrino mass is explained by the celebrated 
seesaw mechanism [0] and the flavor mixing structure is considered to be controlled by 
the Froggatt-Nielsen mechanism 0. There are many works in which the Abelian flavor 
symmetry is discussed ||. On the other hand, there is another experimental suggestion 
on the neutrino oscillation by the Liquid Scintillator Neutrino Detector (LSND) ||. If 
we impose the simultaneous explanation of the result together with the atmospheric and 
solar neutrino deficits, it has been well-known that three different values of the squared 
mass difference are necessary. Then four light neutrinos including a sterile neutrino (y s ) 
are required flOf . Various models of the sterile neutrino can be found in refs. [10-14]. 
Following the recent Super-Kamiokande analysis of the solar neutrino, the explanation 
of the solar neutrino problem based on the v e -v s oscillation seems to be disfavored ||. 
It suggests that the (3+l)-neutrino spectrum might be a more favored scenario for the 



neutrino mass hierarchy than the (2 +2) -scheme [15 



In this paper we consider a neutrino mass matrix in a 3(vl+Vr) framework by using the 
seesaw mechanism. However, being different from the ordinary seesaw models our model 
contains a light right-handed neutrino as a result of the special texture of a right-handed 
Majorana neutrino mass matrix. Although there are the similar works in this direction, 
in most of them it is necessary to introduce the Majorana masses for the left-handed 
neutrinos in order to obtain simultaneously the required values of the mass eigenvalues 
and the flavor mixing angles as it can be found, for example, in O, [TJJ. It means that 



an introduction of a new triplet Higgs field might be necessary. In the present model we 
only need the Dirac neutrino masses and the right-handed Majorana neutrino masses if 
we assume a special but simple texture in both of them at tree level. The model seems 
to have less parameters as compared to the previous ones. 

One of the interesting points of the model is that the large mixing angle MSW solution 
for the solar neutrino problem can be consistently accommodated in the same way as other 
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solutions |16|, [17|]- The LSND result might be also explained if we take an appropriate 



solution for the solar neutrino problem [|I7]]. Moreover, it is interesting that this scenario 
for the neutrino mass matrix could also be embedded into the GUT scheme by introducing 
a suitable flavor symmetry. Such an example in the SU(5) model will be constructed by 
fixing the charge assignment of quarks and leptons for that symmetry. 

The organization of this paper is as follows. In section 2 we define our model and 
discuss its various phenomenological features in the case that the charged lepton mass 
matrix is diagonal. In section 3 we consider the embedding of the scenario into the SU(5) 
GUT scheme. We discuss the realization of the required form of the mass matrix in the 
basis of the Froggatt-Nielsen mechanism. The flavor structure in the quark sector is also 
discussed here. Section 4 is devoted to the summary. 

2. A model of neutrino mass matrix 

We consider a model defined by the following neutrino mass terms which are different 
from the usual seesaw model in the 3(^+7/^) framework: 

- £ mass = H m v<*N v v a + m pi N P N i + \ M p N p N p + h.c., (1) 

a p=2,3 P=2,3 P=2,3 

where v a is an active neutrino [a = e, /x, r) and Np (P = 1 ~ 3) is a charge conjugated 
state of the right-handed neutrino. We make the following assumption for the mass 
parameters in eq. ([[]): 

m 2e = m 2lM = m 2r = fj, m 3e = fj u = m 3r = fj 2 , 

f) ~ f]i ~ fj 2 < 777.21 ~ m 31 ^ M 2 ~ Ms, (2) 

where the mass parameters should be understood as their absolute values, although it is 
not expressed explicitly. A crucial difference from the usual seesaw model is that one of 
the right-handed neutrinos is assumed to be very light and also has very small mixings 
with other heavy right-handed neutrinos. We assume M 23 = in the Majorana mass 
matrix of Np here, for simplicity. Following arguements are not largely changed even 
if we introduce the non-zero M23. Under this assumption we can integrate out heavy 
right-handed neutrinos N p and get the following 4x4 matrix as a result of the seesaw 
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mechanising, 

( A B B D \ 

B C C E 

B C C E 

\D E E F ) 

The matrix elements A ~ F are expressed by the model parameters in (0) as 

^2 ~2 ~2 = = -2 



(3) 



M 2 M 3 ' 



M 
??m 2 i 



M 2 



M 3 ' M 2 



M 3 ' 



I" 71 ,T 5 



7)1X121 mmzt 



rn 



21 



+ 



m 



31 



(4) 



M 2 M 3 M 2 M 3 M 2 M 3 

If we define the diagonalization matrix U of the matrix (0) as m^ iag = U T m u U, we find 
that {/ can be written as 

' cos 9 — sin 9 — sin sin 5 + cos 9 sin 7 ^ 
^sin^ -^cos# ~^p2 "^(cos^si 11 ^ + sin sin 7) 



U 



\ 



^= sin 9 ^= cos $ 75 "75 ( cos 9 sin 5 + sin 6* sin 7) 
— sin 7 — sin 5 1 



(5) 



where | sin 'y | , | sin5| -C 1 is assumed and mixing angles are defined by 

2^2B . Dcos9 + V2E sin 9 —D sin 9 + y/2E cos 9 

sin 7 ~ — , smo " 



tan 29 



A-2C 1 ' F 
The mass eigenvalues of m v are expressed as 



(6) 



mi — A cos 9 + v 7 ^ sin 20 + 2C sin 0, 
m 2 ~ A sin - V2B sin 20 + 2C cos 9, 
m 3 = 0, m 4 = F, 



(7) 



where we neglect the contribution from the fourth low and column of m v to mi )2 taking 
account of the fact such as A > B > ^y- and C > 4r. Here we should note that in 
1 It should be noted that the number of light sterile neutrinos is restricted at most to one in the 
present scenario. We obtain a 3 x 3 matrix if all of Np are heavy. Even in such a case as far as we 
assume a proportional relation between {m\ a ) and (rri2a) as vectors whose components are labeled by a, 
the texture for the active neutrinos is the same as eq. (^). Then it can be applied to the explanation of 
the solar and atmospheric neutrino problems in the same way as the following discussion. It is essentially 
the same as the one discussed in ref. Jl8||, although it is derived in the different context. 
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{a, /3) 


(U) 


-AU m U p% U aj U Pj (=A) 




(I, n) 


(1,2) 


\ sin 2 29 


(A) 


(I, III) 


(1,2) 


\ sin 2 29 


(B) 


(II, III) 


(1,3) 


sin 2 9 


(C) 




(2,3) 


cos 2 9 


(D) 




(1,2) 


-isin 2 29 


(E) 



Table 1. The contributions to each neutrino transition process v a — > vp from each sector (i,j) of the 
mass eigenstates. 

this model the violation of the proportional relation between (m2 Q ) and (m^, a ) as vectors 
is crucial to restrict a number of zero mass eigenvalue into one and to control the mixing 
structure. There is a freedom in a choice of two elements of (m^ a ) which are taken to 
be equal in (0). As far as we consider the case in which the charged lepton mass matrix 
is diagonal, it is not important. But when we consider the different situation, it might 
become crucial to the consideration of the oscillation phenomena. If the charged lepton 
mass matrix is diagonal, the above mixing matrix U is just the flavor mixing matrix 
y(MNS) -^j,^ controls the neutrino oscillation. We assume it in the charged lepton sector 
and also no CP violation in the lepton sector. At this stage we cannot determine to which 
flavor each u a corresponds so that we will use the Roman numerals for the subscript a for 
a while. Next we study the features of the oscillation phenomena in the model in order 
to fix the neutrino flavor. 

The transition probability due to the neutrino oscillation v a — > vp after the flight 
length L is well-known to be written by using the matrix elements of (|5|) as 

(Am 2 \ 
~1# L J ' (8) 

where Am 2 - = |m 2 — m 2 | and the weak eigenstate v a is related to the mass eigenstate 
v% by u a = U a ii)i in the basis that the charged lepton mass matrix is diagonal. In Table 
1 we summalize the contribution to each neutrino transition mode (a, (3) from a sector 
(i,j) of the mass eigenstates. As a phenomenologically interesting case, we consider the 
situation that the mass eigenstates v\ and V2 are almost degenerate and the hierarchy 
(m, 3 Cmi ~ m 2 ) ^ m 4 among the mass eigenvalues is satisfied. This corresponds to a 
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well-known reversed hierarchy scenario for the atmospheric and solar neutrino problems in 
the (3+l)-neutrino spectrum |19]. The absolute value of each mass eigenvalue is smaller 



than the ordinarily discussed scenario because of m 3 = 0. Then every neutrino cannot be 
a hot dark matter candidate. If we apply it to explain the atmospheric and solar neutrino 
data, the squared mass difference should be taken as Jl], |2|] 

2 x 10" 3 eV 2 < Am 2 3 ~ Am 2 3 < 6 x KT 3 eV 2 , (9) 
1(T 10 eV 2 < Am 2 2 < 1.5 x 1(T 4 eV 2 . (10) 

A suitable value of Am 2 2 should be chosen within the above range depending on which 
solution is adopted for the solar neutrino problem. 

By inspecting Table 1 we find that the simultaneous explanation of both deficits of the 
atmospheric neutrino and the solar neutrino is possible if we identify the weak eigenstates 
of neutrinos (e, //, r) with (I, II, III). Under this identification the 3x3 submatrix of (5) 
is recognized as the correctly arranged MNS mixing matrix. If we note that m 3 = 
and Am 2 3 ~ Am| 3 are satisfied, we find that the atmospheric neutrino is explained by 

— > v T obtained as the combination of (C) and (D) in Table 1. This explanation is 
independent of the value of sin#. On the other hand, the solar neutrino is expected to 
be explained by v e — > (A) and also v e — > v T (B). In both processes the amplitude 
A(= —4:J2UaiUf3iUajU/3j) is | sin 2 29. Thus if sin 2 29 ~ 10~ 2 , the small mixing angle 
MSW solution (SMA) is realized fll7f . In the case of sin 2 29 ~ 1, it can give the large 
mixing angle MSW solution (LMA), the low mass MSW solution (LOW) and the vacuum 
oscillation solution (VO) (JT7I] depending on the value of Am 2 2 . The CHOOZ experiment 
| 20[] constrains a component U e % of the MNS mixing matrix pl| . It comes from the fact 
that the amplitude A of the contribution to v e — > v x with the squared mass differences 
Am 2 3 or Am 2 3 always contains ?7 e3 . The model is free from this constraint since [7 e3 = 
is satisfied independently of the value of sin 9. 

In order to see the viability of the scenario in a more quantitative way it is useful 
to estimate numerically what kind of tuning of the primary parameters in ([J) and (|2|) is 
required to realize the suitable value for the oscillation parameters. For the convenience 
we introduce the following parametrization for the three light states: 

" i f\ x i fj 2 i 



Mo VM, VM 



6 



-0.690 

OOOOO 
OOOOOOOO 
OOOOOOOOOO 
OOOOOOOOOOO 

-0.700 - ooooooooooo 
■ ooooooooooo 
ooooooooooo 

°jp 

IFooo 

-0.710 - ooooooooooo 
ooooooooooo 
oooooooooo 
oooooooooo 

O LMA OOOOOOOO 
-0.720 - □ SMA OOOOOOOO 

A LOW OOOOO 

vo 



-0.730 



1.390 1.400 1.410 1.420 1.430 



Fig. 1 The scatter plot of possible solutions for the atmospheric and solar neutrino problems in the 
(ei,e 2 ) plane. 

For simplicity, we assume M2 = M3. Then the overall mass scale is determined by fi 
and the hierarchy among the mass eigenvalues is controlled by E\ and 62- When e\ = \[2 
and €2 = — -^=, two mass eigenvalues mi and m 2 are degenerate. Using the fact, we can 
estimate a typical scale of /i from the condition @ as /i ~ 1.8 x 10~ 2 eV. This value 
corresponds to M2 ~ 5.6 x 10 10 GeV for fj ~ 1 GeV. The deviation from these values of 
ei,2 determines the difference between mi and m-i and also the value of sin 9. In Fig. 1 we 
give a scatter plot of the possible solutions for both the atmospheric and solar neutrino 
problems in the (ej, €2) plane. In this figure solutions for both sign of sin26 l are contained. 
Since we consider the reversed hierarchy here, both possibilities are allowed. From the 
figure we find that the SMA, the LOW and the VO postulate the finer tuning of the 
parameters than the LMA to realize the required values of the squared mass difference 
and sin 2 29. 

In the present model we have a light sterile neutrino. Therefore we may have a chance 
to explain the LSND result, if ~0(1) eV is satisfied. We can check whether all of mass 
eigenvalues and various mixing angles quoted in the above discussion can be consistent 
with the LSND explanation by using the required relations (f|), (^|) and ([/]). If we take 
E\ ~ 1.41 and E2 ~ —0.71 as a typical example in Fig. 1 and also assume M2 = M3 and 
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m2i = rrisi, we obtain 

fh ~ 1.41r}, r? 2 0.71r), sin7 ~ 1.21—, sin 5 ~ 0.21-^-. (12) 

m 31 m 3 i 

In order to see the feature related to the LSND we note that the relevant amplitude ^4lsnd 
can be written by using the unitarity of ^/( MNS ) and the relation 1 7711,2,3! <C | ttt. 4 | as 

a snd = 4(y e r s) ) 2 (^r s) ) 2 - (13) 

Then the amplitude can be written by using eq. (|5|) as 

.Alsnd ~ 2(cos6 l sin7 — sin 9 sin S) 2 (cos 9 sin S + sin 9 sin 7) 2 . (14) 

The LSND data require it to be in the range around 1.2 x 1CT 3 for Aitil SND ~ 1 eV 2 . 
Here we should remind that | sin7|, | sin<5| <C 1 should be satisfied under our assumption 
(0). If we take the large mixing angle solutions for the solar neutrino problem, we obtain 
"4-lsnd ~ 7; (sin 2 7 — sin 2 <5) 2 . Here we impose it to take the above mentioned value, we 
find ?7i3i ~ 5.4r) and then sin 7 ~ 0.23 and sin 5 ~ 0.04 by using eq. (|12"D . Moreover, 777,4 



can take a suitable value for the LSND result such as 7774 ~ -^J 1 ~ 1 eV. On the other 
hand, if we adopt the SMA solution and then cos# ~ 1, we have -4lsnd ~ 2 sin 2 7 sin 2 S. 
If we require it to take the suitable value, we find 77731 ~ 3.2r} and 7774 ~ 0.4 eV which 
is too small for the explanation of the LSND data. Taking account of these analyses, we 
find that the inclusion of the LSND result restricts our model to the large mixing angle 
solutions with respect to the solution for the solar neutrino problem. This feature of the 
model might be favorable if we take seriously the recent Super-Kamiokande analysis of 
the solar neutrino ||. However, even in this case we should comment on an influence on 
the big-bang nucleosynthesis (BBN) due to the oscillation processes u^ T v s in the early 
universe. The required values of sin 7 and sin 5 for the explanation of the LSND data 
induce these processes at a large rate. The BBN bound on v^^ — > v s given in ref. 



cannot be satisfied unless we assume the presence of the large lepton number asymmetry 
at the BBN epock f23|. 

In Table. 1 an only remaining contribution (E) to — > v T cannot imply any evidence 
in the short-baseline experiment even in the case of sin 2 29 ~ 1 since At77 2 2 is too small. 
However, this mode may be relevant to the long-baseline experiment in the case of At77 2 2 ~ 
10~ 4 eV 2 which corresponds to the LMA solution of the solar neutrino deficit. We show 
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Fig. 2 The transition probability Piy^ — > v x i^\) as a function of the flight length L km. We assume 
E =1 GeV, Am 2 13 = 3.5 x 1CT 3 cV 2 and Am 2 12 = 1(T 4 eV 2 . 

the effect of the mode (E) on the Piy^ — *> v x ) in Fig. 2. The dashed line comes from the 
modes (C) and (D) which correspond to the ordinary two flavor oscillation — > v T . The 
thick solid line is the one which is obtained by taking account of the contribution of (E). 
In the thin solid line the contribution of (A) which corresponds to u x = u e is also taken 
into account. This shows that it may be possible to discriminate the model from others 
in the long-baseline experiment such as L > 2000 km. Moreover, the present model may 
be expected to have another experimental signature in the neutrinoless double /3-decay 
24|j . Using eq. ((5|), the effective mass parameter which appears in a formula of the rate 



of neutrinoless double /3-decay can be estimated as 



J2\U ej \ 2 e^m, 



j 



[mi cos 2 6 + m 2 sin 2 9) ~ mi, (15) 



because of the fine degeneracy between rri\ and m 2 . Thus \m ee \ takes the value in the 
range of 0.04 - 0.08 eV which is independent of the value of sin#, that is, the solution of 
the solar neutrino problem. This value seems to be within the reach in the near future 
experiments. 

Before closing this section we order some short comments. In the above analysis we 
assume the reversed mass hierarchy in the (3+1) scheme. As an another hierarchy among 
the mass eigenvalues we can consider the usual one, that is, (m 3 ^ m 2 m{) m 4 ||15|| . 
However, it is found that such a hierarchy cannot be consistent with the experimental 
data in the present scenario if we note that v e cannot be identified with the state / to 
explain the solar neutrino deficit. We should also note that the above result crucially 
depends on the assumption on the charged lepton sector, for which we assumed that its 
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mass matrix was diagonal. As far as the flavor mixing in the charged lepton sector is 
small, the above result seems to be applicable. However, if it is large, the result can be 
largely changed. When we consider the GUT, such situations happen. In the next section 
we study this issue. 

3. Embedding into SU(5) 

We consider an embedding of our neutrino model into the supersymmetric SU(5) 
GUT. In that case the charged lepton mass matrix can be related to the neutrino mass 
matrix through the group theoretical constraint. Thus we cannot assume the small flavor 
mixing in the quark sector independently of the lepton sector. The result obtained in 
the previous section may be modified if we embed our scenario into the GUT scheme. 
In order to control the flavor mixing structure we adopt the Froggatt-Nielsen mechanism 
and introduce Abelian flavor symmetries U(l)j? 1 xU(1)f 2 . The symmetries are assumed 
to be broken by small parameters A and e so that Yukawa couplings inducing the fermion 
masses are suppressed by both powers of A and e. As a result the fermion mass hierarchy 
is produced. We take a model discussed in || as such a typical example and modify it to 
embed our scenario for the neutrino mass into it. 

In the SU(5) GUT quarks and leptons are embedded into the representations of SU(5) 
as follows, 

10 3 (g,w c ,e c ), 5*3(d c ,£), l3u c . (16) 
We assign the charge of U(l)j7 1 xU(l)j7 2 to each representation in the following way f5|: 



10 : 


(3,2,0), 


(0,0,0), 


5* : 


(c,0,0), 


(0,0,0), 


1 


(0,0,0), 





where the numbers in the parentheses represent the charges given to each generation 

and c, a and (3 are non-negative integers. The ordinary doublet Higgs fields H\ and 

i?2 in the minimal supersymmetric standard model are assumed to have no charge of 

11(1)^x11(1)^. In addition to these fields we introduce SU(5) singlet fields S\ and S 2 

which have the charges (—1,0) and (0, —1) of the flavor symmetries, respectively^]. The 

2 We may need to introduce some fields to cancel the chiral anomaly of the flavor symmetry if it is a 
non-anomalous gauge symmetry However, we do not go further into this issue in the present paper. 
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symmetries control the flavor mixing structure by regulating the number of fields S\ and 
5*2 contained in each non-renormalizable term. If the singlet fields Si and 5*2 get the 
vacuum expectation values (Si) and (S2), the above mentioned suppression factors for 



the Yukawa couplings can be realized as the power of A = and e = ffi-. Here M v \ is 



the Planck scale. 

Using the Abelian flavor charges introduced above, we can obtain the quark and lepton 
mass matrices in the following form: 



( A 6 A 5 A 3 
A 5 A 4 A 2 
A 3 A 2 1 



A c e 



(H 2 ), 



M e ~ 



( A 3+c A 2+c A c \ 
A 3 A 2 1 
A 3 A 2 1 

I A 3+c A 3 A 3 \ 

A 2+c A 2 A 2 
A c 1 1 



(Hi), 



V 



(Hi), 



( e 2a 



M, 



e a+f3 e a+f3 \ 



\ 



c2d 



,2/3 



M, 



J 



where M is the mass scale relevant to the origin of the right-handed Majorana neutrino 
mass. Dirac mass matrices are written in the basis of ipRrriDipL- We do not consider the 



CP phases here. In the mass matrices fll8| ) we abbreviate the order one coupling constants 
by using the similarity symbol. We should note that M v and Mr in QT£| ) can have a similar 
texture to the one defined by eqs. (0) and (0) up to the implicit coefficients of order one 
as far as a > is satisfied^. At least in the case of c = and 1 which is assumed in 
the following discussion, we can make M u satisfy the condition @ by tuning the order 
one coefficients. Thus we can have the similar mass matrix to eq. (H) as a result of the 
seesaw mechanism, although there is a non-zero element M23 in Mr differently from the 
one defined by eq. ([!]). Its diagonalization matrix U can be considered to have the similar 
form as eq. (|5]). Their difference comes only from the definition of the matrix elements 
A ~ F as we will see it below. 

In the quark sector the mass eigenvalues and the CKM matrix elements can be found 
In this context the order one coefficients assumed here may be allowed to be considered in the range 
(v^t;) if e< A for M v . 
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after some inspection as 

m u : m c : m t = A 6 : A 4 : 1, : m s : = A 3+c : A 2 : 1, (19) 

Ks~A, V ub ~\\ K fe ~A 2 . (20) 

On the charged lepton sector we can know the mass eigenvalues by noting the SU(5) 
relation such as Mj = The ratio of mass eigenvalues is the same as the one of the 
down quark sector and then 

m e : m M : m T = A 3+c : A 2 : 1. (21) 

The result has some different features from the ones presented in ref. || in the down quark 
and charged lepton sectors. It comes from the charge assignment for 5* which is needed to 
realize the Dirac neutrino masses defined by eqs. (P and (0). If we assume A ~ 0.22, these 
results seem to describe the experimental data in a qualitatively favorable way, except 
for m e and m u which are predicted to be too large, in particular, in the case of c = 0. 
This is the common fault known in the scheme based on the Abelian flavor symmetry and 
its similar charge assignment to the one given in fllTD - We cannot overcome it without 
something new. 

We define the diagonalization matrix U of the charged lepton mass matrix in a basis 
that U^M\M e U is diagonal. Then U can be approximately written as 



U 



^cos£-^sin£ 
-72 COS £-76 sin £ 



75 sin £+76 COS £ 
-72 sin £ + 76 COS ^ 
76 COS ^ 



i 



o\ 



u 



V 

( cos£ 
-75 sin £ 











1 





1 











(22) 



V 



o 



sin£ \ 
^cos£ 

1 J 



1 

V2 



(23) 



The hierarchical structure ( PT| ) of the mass eigenvalues requires a mixing angle £ to be 
sin£ ~ 0(\). In the neutrino sector we need to determine the finer structure of the 
Dirac neutrino mass matrix to be suitable for the charged lepton mass matrix given in 
(111) from a viewpoint of the explanation of various data for the neutrino oscillations. 
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For that purpose we should remind that there is a freedom in the choice of two elements 
of Dirac neutrino masses (m 3a ), which are taken to be equal by tuning of the order 
one coefficients. After some investigation we find that it seems to be favorable to take 
m 3e = m 3fl = fji, m 3r = fj 2 instead of the one given in eq. (Q) 4 . Under this assumption 
the mass matrix of the light neutrinos can be written as 

( A A B d\ 

A A B D 

B B C E 

\D D E F ) 

and the matrix elements A ~ F are defined by 

77771 



m. 



(24) 



A = jf_ rg_ 
M 3 M 2 



V , V1V2 f)(Vi + V2) 



'M23' 



B = -!— + 

M 3 M 2 



M 



23 



c 

E 



M 3 M 2 M 23 ' 



D 



f)m 21 f|im 31 f)m 3l + r] 1 m 2 i 



M 3 



Mo 



M< 



2?, 



77777,21 r] 2 m 3l fjm 3l + r] 2 m 21 



rn 



21 



+ 



rn 



31 



m 21 m 31 



(25) 



M 3 M 2 M 23 M 3 M 2 M 23 

where M~ x = M a j{M 2 M 3 — M 23 ) and we partially use the notation in eq. (||). The mass 
eigenvalues of are 

mi ~ 2A cos 9 + V2B sin 29 + C sin 6, 
m 3 ~ 2 A sin 2 9 - y/2B sin 29 + C cos 2 9, 

m 2 = 0, m 4 = F, (26) 

where we again neglect the additional contributions to 777^3 because of the same reason as 
the one in the previous section. The diagonalization matrix U is rearranged from eq. @ 
because of the change in the choice of fji and f] 2 . Using the modified U and eqs. ([22]) and 
(p3[), the MNS matrix of the lepton mixing defined by \/( MNS ) = JJ T U is calculated for 
both values of c as 



V, 



(MNS) 



c=0 



sin£ 
cos £ 

(2) 



- sin 7 



4 It means that fj ~ 771 ~ A c e /3 ~ e 13 and 
one coefBcients. 



cos £ 
sin£ 







sin£ a\ 



■> + 

V3 
f (2) 



\ 



- sin 5 



ff3 
1 



(27) 



iP . The difference among them comes from the order 
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ta( mn S) „ 
c=l 



2 2~~ ffl 2 

^sin£ + ^cos£ ±cos£ sin£ + cos£ a 3 



(28) 



^ — sin 7 — sin 5 1 y 

where we use definitions 

aj = f n sin 7 + t> i3 sin 5, = cos 8 ± \/2 sin 8, f±^ = \[2 cos 8 ± sin 6 1 , (29) 

and represents the r/'-element of the corresponding \A MNS ). To derive these expressions 
we use I sin7|, | sin<5| <C 1 and neglect higher order terms of them. The mixing angles 
8, 7 and S in this case are defined as 

2V2B . y/2D cos 8 + E sin 8 . . -V2Dsm8 + E cos 8 , . 
tan 20 = 2A _ C , sm 7 ~ , sm5 ~ - . (30) 

Now we study the oscillation phenomena in both cases in more detail. First we consider 
the case of c = 0. Taking account that m 2 = and l / r2 4NS ' ) = 0, the reversed hierarchy 
scenario cannot be adopted from a viewpoint of the atmospheric neutrino problem. We 
must assume the normal hierarchy (|m 2 | ^ \mi\ <C \m-i\) <C \m^\ in the (3+l)-scheme in 
order to realize Am^ ~ Am s 2 olar and Am^ ~ Am 2 3 — Am 2 tm . Using eq. (|27D , we find 
that the amplitude for — > v T is 

A = _4V;( 1 MNS V r ( 1 MNS V A ;3 MNS V T ( 3 MNS) = -(V2 - - tan 2df cos 2 28 cos 2 £. (31) 

After some investigation we find that it suggests that tan 28 < and cos 2 8 ~ 1 should be 
satisfied for the explanation of the atmospheric neutrino problem. We should also remind 
the fact that sin£ ~ 0(A). Although the large value of |sin#| such as 0.95 can satisfy 
the bound from the atmospheric neutrino, it seems to be disfavored by the solar neutrino 
data. Thus the atmospheric and solar neutrino problems can be explained by — > v T 
corresponding to Am 2 3 and v e — > corresponding to Amf 2 , respectively. The situation 
for the solar neutrino problem is different from the case in the previous section. Only 
the SMA solution is allowed in the present case since K ( 2 MNS) = and m 2 = make the 
contribution of v e — > v T to solar neutrino deficit zero. It originally comes from the non- 
diagonal structure of the charged lepton mixing matrix fl2"2"|). Since the present neutrino 
mass matrix induces the large mixing between v e and by itself, we need the small 
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mixing in a corresponding place of the charge lepton sector to realize the large mixing 
solution for the solar neutrino. However, it is not satisfied there in this case. 

In order to see the viability of the model quantitatively we need to check numerically 
the consistent realization of both data of the atmospheric and solar neutrino observa- 
tions. In this study, for simplicity, we assume M 2 3 = and M 2 = M 3 here, although 
only the symmetries 11(1)^x11(1)^ cannot verify the former one. Then we can use the 
parametrization (|TT|). In Fig. 3 we give the scatter plot of solutions for both the atmo- 
spheric and solar neutrino problems in the (61,62) plane assuming —0.42 < tan 28 < 0. 
Using the figure, we can find the typical values of the primary parameters in the model 
by using eqs. (|25|) , ( |26|) and (|30|) . As an example, if we take e\ ~ 3.0 and e 2 ~ —0.45 from 
Fig. 3, we can obtain 

m~3M, fj 2 0.45?}, tan2#~-0.05, sin 7 ~ 2.8—, sin 5 ~ 0.35-^-. (32) 

TO31 m 31 

If we assume cos£ = 0.98 and sin£ = 0.2, in this case the MNS matrix becomes 

/ -0.12 -0.98 0.16 -0.12 SU17 + 0.16 sin <5 \ 



(MNS) 



0.59 -0.20 -0.78 0.59sin7-0.78sin<5 
0.80 0.60 0.80sin7 + 0.60sin5 
- sin 7 — sin S 1 J 



(33) 



The CHOOZ constraint on V e 3 MNS ' ) is satisfied. The LSND result may have a chance 
to be again explained because of the exsistence of one light sterile neutrino. In order 
to see it we study the relevant amplitude -4lsnd which is estimated by eq. (|H|). If we 
require ^4lsnd ~ 1-2 x 10~ 3 , we obtain 777,31 ~ 4.7?7 and then WI4 ~ 0.1 eV where we take 
/i ~ 2.6 x 10~ 3 eV. This is too small to explain the LSND data. Although the larger 
value of ^4lsnd induces the smaller value of m 4 , we cannot find a favorable result for the 
LSND within the present freedom. In the present case [i and V^i MNS ' ) tends to take small 
values by the requirement of the atmospheric and solar neutrinos. As a result of this 
general feature, 777,4 takes a small value compared to the required value by the LSND. 
The effective mass parameter \m ee \ for the neutrinoless double /3-decay can be estimated 
as \m ee \ ~ f/g 4 m 4 ~ 0.16/x It is too small as compared with the value expected to be 
reached by the near future experiment. 

It is also useful to note that the above values of the primary parameters of the model 
are realized through the suitable charge assignment of a and j3. In order to show such an 
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Fig. 3 The scatter plot of possible solutions for both of the solar and atmospheric neutrino problems in 
the (ei, £2) plane. Requiring —0.42 < tan 28 < 0, the scale parameter fi is taken as 7.7 x 10~ 3 eV for the 
LMA and 2.6 x 1CT 3 cV for the SMA. 

example we take (H 2 ) ~ 100 GeV, M~4x 10 15 GeV and e ~ 10~ 2 . Then if we assign 
a = 6 and (3 = OQ, we can check that all required quantities except for .4lsnd are realized 
in the suitable range discussed numerically above up to the order one factors. In this case 
we have 777,4 ~ 1 eV. 

Next we treat the case of c = 1. Also in this case the reversed hierarchy cannot induce 
the sufficiently large amplitude A for — > v T using the modes with Am 2 2 and Am| 3 . We 
need adopt the normal mass hierarchy also in the present case. The relevant amplitude 
for — > v T is estimated as 

(34) 

where we take account of sin £ 1 in the estimation. If cos 2 9 ~ 1 is satisfied, the above 
amplitude can be suitable to the atmospheric neutrino problem. The contribution to the 
solar neutrino deficits comes from u e — > and v e — > v T with Am 2 2 . Their combined 
amplitude is almost equal to one and the large mixing angle solution is realized. The 
value of cos 2 9 is fixed to generate the large mixing for the explanation of the atmospheric 
5 The more general condition is e a = 1CP 12 . Since e can be taken as a very small value, it can be 
consistent even if the coefficients which are considered to be of order one are rather large. 
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neutrino and it also results in the large mixing between v e and v^. On the other hand, 
in the charged lepton sector c = 1 makes the mixing between e and /i small so that we 
can have a large mixing angle solution for the solar neutrino. Also in this case we give 
the scatter plot of the LMA solutions in Fig. 3 by assuming M 23 = and M 2 = M 3 and 
using the parametrization ([11]). The LOW and VO solutions seem to be difficult to be 
realized since the required Amf 2 should be much smaller than the LMA. Using Fig. 3 we 
can again find the typical values of the primary parameters in the model. As an example, 
if we take e\ ~ 1.75 and e 2 ~ —0.75, we can obtain 

77i~1.75r), r/2 0.7597, tan2# 0.13, sin 7 ~ 1.9—, sin 5 ~ 0.26—. (35) 

777,31 rn 31 

These fix the MNS matrix in the present case as 



T/ (MNS) 
V c=l 



t 


0.62 


-0.79 


-0.10 


0.62sin7-0.10sin<5) 


\ 




0.55 


0.50 


-0.67 


0.55 sin 7 — 0.67 sin S 






0.58 


0.35 


0.73 


0.58 SU17 + 0.73 sin 5 




V 


— sin 7 





— sin S 


1 


/ 



(36) 



We can see that the CHOOZ constraint on V'J 3 MNS ' ) is satisfied in (|36|) . In order to see the 
possibility to explain the LSND result we impose v4lsnd ~ 1-2 x 10~ 3 . Then by using 
eq. (|I3|) and fi ~ 7.7 x 10~ 3 eV, we obtain m 31 ~ 7.8?} and then m 4 ~ 0.93 eV which seems 
to be in the suitable region for the LSND. We also have sin 7 ~ 0.25 and sin S ~ 0.03 
which are consistent with our assumption for sin 7 and sin 5. These analyses show that 
this case can explain all neutrino oscillation data including the LSND. It is also interesting 
that the effective mass \m ee \ for the neutrinoless double /3-decay can have rather large 
value because of the 7774 contribution. In fact, using the above numerical values it can 
be estimated as |77? ee | ~ |?7 e 4| 2 77?4 ~ 0.02 eV, which may be a promising value from the 
experimental viewpoint. The value of the charges a and (3 adopted for the c = case 
is also applicable to the present case to realize the model parameters in a favored region 
up to the order one coefficients as far as we take the same values of {H 2 ) ~ 100 GeV, 
M ~ 10 15 GeV and e ~ 10~ 2 . However, the present case needs more complicated structure 
of the order one coefficients as compared with the c = case. 

Finally, we should note that in our scheme the required equality among (m 2a ) and 
also among (777 3a ) in eq. (0) might not be necessary to be satisfied exactly. The allowed 
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deviation should be quantitatively investigated since it is related to the estimation of the 
magnitude of order one coefficients. 

4. Summary 

We have proposed the scenario for the neutrino mass and mixing based on the seesaw 
mechanism in the 3(u L + u R ) framework. By assuming the special texture for the right- 
handed Majorana neutrino mass matrix and the Dirac mass matrix we could obtain a 
model with four light neutrino states including a sterile neutrino. One active neutrino 
is massless and others can have the masses which are suitable for the explanation of 
the atmospheric and solar neutrino deficits, and also the LSND result. We studied two 
different cases specified by the diagonal charged lepton mass matrix and the non-diagonal 
one which is obtained by embedding of our neutrino mass matrix into the SU(5) GUT 
scheme. 

In the former case the so-called reversed mass hierarchy scenario has been adopted. 
Every known solution for the solar neutrino problem could be realized by tuning the Dirac 
mass matrix of neutrinos. It is interesting that the large mixing angle MSW solution 
can be most easily realized as compared to other solutions. Moreover, if we impose the 
explanation of the LSND on it, the solution for the solar neutrino problem is restricted to 
the ones with the large mixing angle. The difference from the two flavor oscillation could 
be expected to be observed in the — > v T using the long-baseline experiment with the 
flight length more than 2000km. The neutrinoless double /5-decay might also be accessible 
if the experimental bound is improved to the level of \m ee \ ~ 0.04 - 0.08 eV. 

In the latter case the Froggatt-Nielsen mechanism has been applied to control the flavor 
mixing. We introduce the Abelian flavor symmetries U^)^ xU(1)_f 2 whose factor groups 
are assumed to have the different breaking scale. The non-trivial charge assignment of 
U(1)f 1 is used only for the 10 and 5* fields of SU(5) and the right-handed neutrino 1 is 
assumed to have only the charge of U(l)i? 2 . Under this setting we studied the features 
of the mass and the mixing in both quark and lepton sectors for the two types of the 
charge assignment. We found that it could generate the mass eigenvalues and the flavor 
mixings for the quark sector in a qualitatively satisfactory way. If we give the suitable 
flavor charge to the right-handed neutrinos, our neutrino scenario can be also embedded 
into the SU(5) scheme consistently. Although the reversed mass hierarchy is disfavored in 
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both charge assignments, the ordinary mass hierarchy presents a consistent explanation 
of all data of the known neutrino oscillation observations. For the solar neutrino problem 
only the SMA solution or the LMA solution is allowed in each case. However, if we impose 
the explanation of the LSND result, only the LMA seems to be favored. 

In this paper we assumed that the tuning of order one coefficients could always be 
allowed. Although in our scenario the mild tuning of order one coefficients is very crucial, 
we cannot say anything on its origin at the present stage. 
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